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ABSTRACT 

Integral field spectroscopy of 11 type-Ib/c supernova explosion sites in nearby galaxies has been 
obtained using UH88/SNIFS and Gemini-N/GMOS. The use of integral field spectroscopy enables us 
to obtain both spatial and spectral information of the explosion site, allowing the identification of the 
parent stellar population of the supernova progenitor star. The spectrum of the parent population 
provides metallicity determination via strong-line method and age estimation obtained via comparison 
with simple stellar population (SSP) models. We adopt this information as the metallicity and age of 
the supernova progenitor, under the assumption that it was coeval with the parent stellar population. 
The age of the star corresponds to its lifetime, which in turn gives the estimate of its initial mass. 
With this method we were able to determine both the metallicity and initial (ZAMS) mass of the 
progenitor stars of the type lb and Ic supernovae. We found that on average SN Ic explosion sites are 
more metal-rich and younger than SN lb sites. The initial mass of the progenitors derived from parent 
stellar population age suggests that SN Ic have more massive progenitors than SN lb. In addition, 
we also found indication that some of our SN progenitors are less massive than ~ 25 M , indicating 
that they may have been stars in a close binary system that have lost their outer envelope via binary 
interactions to produce Ib/c supernovae, instead of single Wolf-Rayet stars. These findings support 
the current suggestions that both binary and single progenitor channels are in effect in producing type 
Ib/c supernovae. This work also demonstrates the power of integral field spectroscopy in investigating 
supernova environments and active star forming regions. 
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Subject headings: supernovae: general — stars: massive 



1. INTRODUCTION 

A great number of supernovae (SNe) have been ob- 
served and new ones are being discovered owing to ongo- 
ing transient surveys and monitor ings. Despite the large, 
growing number of observed events the real nature of the 
progenitor of each type of supernova is still not very well 
understood. 

Recently, a number of core-collapse supernovae (CC- 
SNe) have had their progenitor st ar detected in h i gh res - 
olution pre-explosi o n images (e.g. [ Van Dvk et al. (2 0121) : 
iFraser et all (|2011f ): lElias-Rosa et al.l ([201 see lSmarttl 
(|2009[ ) for a review). The detection of the progenitor 
in several passbands allows the comparison between ob- 
served colors and luminosity of the star with stellar mod- 
els. While direct detection is obviously desirable, in most 
cases only upper luminosity limits could be derived due 
to non-detection. This, however, is still useful in con- 
straining the progenitor properties. Comparison with 
theoretical models yields the evolutionary stage and po- 
sition in the Hertzprung-Russel diagram of the star, from 
which the in itial mass cou l d ultim ately be derived. Using 
this method ISmartt et al.l (|2009f ) concludes that SNe II-P 
arise from red supergiant stars with initial mass between 
8.5±1.5 and 16.5 ± 1.5 M . 

While the progenitors of SNe II-P have been convinc- 
ingly detected and characterized, currently no Ib/c pro- 
genitor has been detected. It is widely believed that since 
SNe Ib/c are deficient in hydrogen and helium, the most 
likely candidate for their prog enitors would n aturally be 
the Wolf-Rayet (WR) stars (jCrowtherl [2QQ7h . Numer- 
ous stellar evolution models have been proposed to char- 
acterize the progenitor stars of different S N types (e.g. 
Hege r et al.l 120031 : lEldridge fc Toutl 120041 : iGeorgy et all 
2009) but these predictions are still not very well tested 
with observational data. Progenitor discoveries could 
provide useful comparisons but the number of data points 
is still too small and not very reliable, considering a num- 
ber of uncertainties such as metallicity assumption and 
also the sensitivity of the method to errors in luminos- 
ity (determined using host galaxy distance) and the end 
state of the evolution of the exploding star. Furthermore, 
progenitor discoveries ideally should be followed-up by 
post-explosion observation to ensure that the purported 
progenitor has indeed gone. However, this has not been 
done for the majority of the progenitors thus future con- 
firmation of their disappearance is still needed. 

Initial mass of a progenitor star is generally considered 
as the most influential parameter driving the evolution 
of the star. Hydrogen-deficient (type-lb, Ic) SNe require 
progenitors with high mass-loss rate to remove the outer 
hydrogen e nvelope of the sta r, in which metallicity- driven 
wind (e.g. IVink et all l2001f ) or other mechanism could 
play a significant role. Therefore, it is also important 
to investigate the metallicity of the progenitor star in 
addition to the initial mass. This dependence on metal- 
licity affects the number rat io N Ib / c /Nu as sug gested 
by several previous studies (jPrantzos fc BoissTerl |2003; 
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Boissier & Prantzos 2009; Eldridge et al. 2008). Another 
possible mech anism of mass loss is via close binary inter- 
action (e.g. Podsiadlowski et al. 1992). Extending the 
work of Smartt et al. (2009) for SN Ib/c progenitors, 
lEldridge et al.l {2013) suggest that the non-detection of 
SN Ib/c progenitors may be explained by the idea that 
probably most of SN Ib/c progenitors arise from in- 
teracting binaries with mass < 20 M . An example 
where binary scenario provides a good alternative expla- 
nation to the massiv e single star progenitor is given by 
iCrockett et all (j2007D in the case of SN Ic 2002ap. From 
analysis of pre-explosion images, they suggest that the 
progenitor of SN 2002ap was probably initially a 30-40 
M star with twice the standard mass-loss rate, or in a 
b inary system as a 15- 20 M star. 

I Anderson et al.l (j2Q10h showed that SNe Ib/c occur in 
higher metallicity environments than SNe II. High metal- 
licity massive stars tend to explode as SNe Ib/c, thus 
higher number ratio N Ib / c /Nu is expected in higher- 
met allicity environment. While I Anderson et al.l {2010) 
suggest that SN lb and Ic have similar metallicities, 
iModjaz et al.l (j2011f ) argue t hat SN Ic have h i gher m etal- 
licity th an lb. The results of ILeloudas et~aT1 (j2011l ). and 
recently iSanders et al.l (|2012f ) , also show indications for 
higher metallicity for Ic, but those were not statistically 
significant. In contrast with this work, all those previous 
studies mentioned above used slit spectroscopy to obtain 
spectra of the SN explosion sites, from which the site 
metallicities were derived. 

Considering the limitations of finding SN progenitors 
in high-quality pre-explosion images, several alternative 
methods have been devised to obtain insights into the 
nature and possibl e mass range of the prog enitors. A 
statistical study by lAnderson fc Jamesl (j2008[ ) using cor- 
relations of SN positions within host galaxies with Ha 
emission showed that the association of SNe with recent 
star formation follows a SN Ic > lb > II sequence, in- 
dicating a grouping of decre asing progenitor mass. This 
res ult is furth e r rein forced in lAnderson et all {2012); but 
see lCrowthe"rl (j2013l ) for a discussion on how SN associ- 
ation with a nearby HI I only con strains the progenitor 
mass weakly. ILeloudas et al.l (j2011f ) suggest that SN Ib/c 
progenitors are more massive than 25-30 M from the 
age of the youngest stellar population at the explosion 
site, estimated using Ha emission line equivalent width. 
However, in that case it is likely that multiple stellar 
populations have been observed and treated as a single 
one since the study used integrated light over a large 
region within the spectrograph slit, in some cases ~ 1 
kpc. Some of their targets are also not the exact explo- 
sion sites but other regions several kpc away in the host 
galaxy. Therefore, consequently the real parent popula- 
tion might has been contaminated by other populations 
of possibly different ages and metallicities. This problem 
i s als o prese nt in similar s tudies such as IModjaz et al.l 
(j2011f ) and ISanders et al.l (j2012[ ). Here we endeavour 
to refine the results of those efforts by resolved spec- 
troscopy of the explosion site, with typical spatial resolu- 
tion of around 70 pc ( 1" at 15 Mpc). A simi lar approach 
has been employed bv lGogarten et all (|2009l ) to estimate 
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the initial mass of the progenitor of the optical transient 
NGC 300 OT2008-1. They probed region -50 pc around 
the transient, and determined the initial progenitor mass 
of 12-25 M from the analysis of the stellar population 
at the site. 

In this paper we report the results of our study of 
11 nearby Ib/c SN explosion sites, in which the par- 
ent stellar populations were detected and used to de- 
rive progenitor star's mass and metallicity. We will use 
the terms star cluster, HI I region, and OB association 
interchangeably in this paper to refer to t he stellar pop- 
ulatio n. As s tars are born i n clus ters (jLada fc Ladal 
2003, also see iBressert et al.l (j2010h who suggest that 
this depends on the adopted cluster definition, which 
may change the fraction of stars born in clusters to 
be ~ 45-90%) and there is evidence that all mas- 
sive stars may be bor n within clustered environments 
(jPortegies Zwart et al.ll2010h . it is reasonable to derive 
a coeval star's age and metallicity from its parent stel- 
lar population. It is not expected that the progenitor 
star would have wandered far from the birthplace clus- 
ter since they would not have much time to travel before 
exploding as a supernova - thus the parent cluster may 
be studied to extract useful information on the coeval 
SN progenitor star. Typical velocity dispersion inside 
a star cluster is a few km/ s (jBastian fc Goodwin! 2006; 
iPortegies Zwart et al.|[20To[ ) , which corresponds to about 
a few pc/Myr. Therefore a SN progenitor which belongs 
to such cluster, if unbound, would only travel typically 
a few tens of par sec during its several Myr lifetime. By 
studying the host cluster we were able to derive initial 
mass and metallicity of the SN progenitors. 

We organize the paper as follows. Observation, data 
reductions, and the subsequent analysis is presented in 
Section [2] Descriptions of each SN explosion sites are 
given in Section [3l We discuss our findings in Section |4] 
and finally a summary is presented in Section [5] 

2. DATA ACQUISITION AND ANALYSIS METHOD 

We selec ted our samples fro m the Asiago Supernova 
Database (jBarbon et al.l fl999h . SNe in host galaxies 
with radial velocity of 3000 km/s or less were selected 
and subsequently careful visual inspections were per- 
formed to DSS and SDSS images of the explosion sites 
using ALADIIsQ. We use equatorial coordinates from the 
Asiago database to locate the SN positions in the im- 
ages. In addition, we also used the publishe d images 
of the SN environment/light echo study by Boffi et al. 
(1999) to compare with and identify some SN host clus- 
ters. SN classification follows the Asiago database clas- 
sification; in this paper we present the study of our 
S Ne Ib/c samples whi l e the type-II SNe are presented 
in iKuncaravakti et al.l (|2013l Paper II). 

SNe showing close association with a bright knot were 
selected as observation targets. We interpret these knots 
as the parent stellar population of the SN progenitor star. 
The utilization of broadband images ensure that there is 
no selection preference towards very young stellar pop- 
ulations, as opposed to if we carried out our selection 
using Ha or UV images. Therefore, we expect very little 
age bias to be present in our sample explosion site. 

1 http : //aladin . u- strasbg . f r/aladin . gml 



We observed 7 of our targets in 2010-2011 using 
the Super Nova Inte gr al Field Spectrog raph (SNIFS; 
lAldering etldl (|2002D . lLantz et all (|2QQ4f )) mounted at 
University of Hawaii 2.2 m telescope (UH88) at Mauna 
Kea. Table Q] contains the observations of our SN site 
targets. The positional uncertainty of each SN is also 
showed; the reasons for each estimates are given in the 
description of each explosion sites. 

The integral field spectrograph SNIFS employs a fully- 
filled 15 x 15 lenslet array in the integral field unit (IFU), 
covering 6.4" x 6. 4" wide field of view (FoV). The cor- 
responding spatial resolution is 0.43 arcsec per spatial 
pixel (spaxel). SNIFS has a photometric channel that 
is used to accurately place objects in the IFU. The po- 
sitional accuracy with which a given (RA, Dec) coordi- 
nates is placed on the IFU is around 0.2 arcsec. The 
scale and rotation of the SNIFS IFU is well measured 
thus should not produce much additional astrometric un- 
certainty. The spectrograph consists of two arms oper- 
ating simultaneously to cover the whole optical wave- 
length range at R ~ 1000: the blue channel covers 
3300 - 5200 A and the red channel covers 5100 - 9700 
A. SNIFS is operated remotely with a fully-dedicated 
data reduction pipeline. The final pipeline result is 
wavelength- and flux-calibrated (x, y. A) datacubes. The 
descripti on of the pipe l ine is similar to what is presented 
in §4 of iBacon et al.l (j20Qlh and is outlined briefly in 
lAldering et al.l (|2006f ). 

Observation of the other 4 targets was carried out 
using the Gem ini Mult i- Object Spectrograph (GMOS) 
in IF U mode (jAllingt on- Smith et al.l 120021 : iHook et al.l 
2004), with the 8.1 m Gemini North telescope at Mauna 
Kea (Program ID GN-2011B-C-6, via Subaru-Gemini 
Time Exchange Program). GMOS was used in one-slit 
mode, providing a 5" x 3.5" field of view sampled by 
the hexagonal lenslet array with one lenslet diameter of 
0.2". We used the B600 grating centered at 5400 A, and 
also took identical exposures with the grating centered 
at 5450 A in order to perform spectral dithering to elim- 
inate cosmic rays and cover the gaps of the three CCDs 
utilized by GMOS. This enables us to observe continu- 
ous spectra from 4000 A to 6800 A at R ~ 1700, thus 
covering and resolving prominent spectral lines such as 
Ha, H/3, H 7 , [OIII]AA4959,5007, [NII]AA6548,6583, 
and [Sll]AA6716,6731. Subsequent data reduction was 
carried out using the Gemini package in IRAF0. The 
data was reduced to remove instrumental signatures and 
then calibrated in wavelength and flux to produce final 
(x, A) datacubes. 

Once the raw SNIFS and GMOS data have been fully 
reduced into (x, y. A) datacubes, object spectrum extrac- 
tion and analysis were done using IRAF. The datacubes 
were treated as stacks of images taken at different wave- 
lengths. The flux density of an object at each wavelength 
'slice' was extracted by performing aperture photometry 
(IRAF/ app hot) on the object at each slice, and arrang- 
ing the values in the wavelength direction to produce a 
spectrum. The seeing FWHM was used as the radius of 
the circular aperture, or in some cases smaller apertures 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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were used due to object proximity to the edge of the 
field of view. Sky in SNIFS data was estimated using an 
annular aperture around the object and subtracted out. 
Usually the annulus is larger than the field width thus 
effectively only a part of it was being used to measure the 
sky. GMOS is equipped by another IFU dedicated for si- 
multaneous sky observation during the science exposures, 
and subtraction was done in the reduction process. 

The positions of the objects in each SNIFS wave- 
length slice were traced, thus compensating the differ- 
entia l atmospheric refraction effect (DAR; Filip penkol 
1982). This demonstrates the advantage of using inte- 
gral field spectroscopy (IFS) compared to conventional 
slit spectroscopy which is strongly affected by DAR, and 
with IFS the effects of DAR could be corrected a pos- 
teriori as has been shown by lArribas et al.l (|1999l ). The 
GMOS dataset was corrected for DAR in the datacube 
stage by shifting each wavelength slices, and aperture 
photometry of individual slices was done with fixed aper- 
ture position for each object. 

The final resulting spectra were then analysed using 
IRAF / splot. Nebular emission lines were measured by 
fitting a Gaussian profile. To measure metallicity, only 
the nebular emission lines are needed thus the stellar con- 
tinuum was removed prior to line measurement by fitting 
a polynomial function. Before measuring line equiva- 
lent width, the continuum was normalized also by fit- 
ting with a polynomial. Metallicity in terms of oxygen 
abundan ce was derived u s ing th e 03N2 and N2 calibra- 
tions bv iPettini fc Pagell (|2QQ4 hereafter PP04). This 
method of metallicity determination uses intensity ratios 
of lines closely spaced in wavelength, therefore not sen- 
sitive against errors introduced by reddening and flux 
calibration. 

Following PP04, we adopt 12 + log ( O/H) = 8.66 as 
the v alue of the solar oxygen abundance (Asol und et al.l 
l200l . In all of the observed fields Ha and [NIIJA6583 
were detected thus enabling N2 determination, while the 
H/3 and [OIIIJA5007 detection required for 03N2 is less 
frequent. In cases where both 03N2 and N2 were de- 
tected we adopt the mean value as the metallicity, other- 
wise the N2 value is adopted. The errors quoted in metal- 
licity in Z unit are the bounds for the highest and lowest 
values of 12+log(0/H) determined from 03N2 and N2. 
In cases where only N2 is available, we assign the error of 
12+log(0/H) as ±0.18 dex, which is the la uncertainty 
in the PP04's N2 calibration determination. These were 
then converted into Z metallicities to give the highest 
and lowest value bounds of metallicity. 

Ages of the stellar populations were determined by 
comparing the observed Ha equivalent width (EW) 
with the theoretical values provided by simple stel- 
lar po pulation (SSP) model Starburst99 (jLeitherer et al.l 
1999) for each appropriate metallicity. We assume an 
instantaneous-burst stellar population with continuous, 
standard Salpeter IMF (a = 2.35). The evolution of 
HaEW with SSP age is presented in Figure [U HaEW 
of a stellar population with continuous star formation is 
also shown for comparison. The evolution of Ha is sen- 
sitive to the star formation history, while IMF variations 
give less important effect especially at SSP age older than 
~ 4 Myr. The error of HaEW were estimated from the 
signal-to-noise ratio of the spectral continuum around 
Ha. We expect the effect of dust absorption in HaEW 
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Figure 1. Time evolution of Ha equivalent width, Starburst99 
SSP. Single-burst model is represented with black line color and 
continuous star formation with purple; both use standard Salpeter 
IMF with a = 2.35, M U p = 100 M . Dotted lines represent single- 
burst solar-metallicity models with different IMFs: red is for a = 
3.30, M U p = 100 M , orange is for a = 2.35, M U p = 30 M . The 
lifetimes of single stars of different initial masses at solar metallicity 
according to Padova models are indicated with vertical grey lines. 

determination to be minimal due to the very small wave- 
length range used for EW measurement. The SSP age de- 
termined from HaEW is adopted as the lifetime of the co- 
eval progenitor star, and comparison with Padova stellar 
evolution models gives the estimated initia l mass of a star 
with such lifetime. We used models from Bressan et al. 
(1993) for solar metallicity (Z = 0.02), and Fag otto et al.l 
(|1994f ) for 0.4 solar metallicity (Z = 0.008). We defined 
observed oxygen abundance of 0.7 (O/H)©, correspond- 
ing to 12+log(0/H) = 8.50, as the dividing line between 
using the solar-metallicity model or the subsolar one to 
determine stellar initial mass. These models extend from 
the zero-age main sequence (ZAMS) to the onset of cen- 
tral carbon burning. This is considered representative of 
the stellar lifetime since the time elapsed from carbon 
burning to SN is only about a few hundred years, neg- 
ligible compared to the ~ 10 6 — 10 T year lifetime of the 
star. 

We compare Starburst99 model of HaEW a gainst 
GALEV (lAnders fc Fritze-v. Alvenslebenl 120031 : 
iKotulla et al.lH2009[ h to check for consistency and 
possible systematic effects. Both SSP models consider 
the contribution of stellar and nebular continuum to 
the output light, but GALEV includes nebular line 
emissions. As for the HaEW value as a function of SSP 
age, Starburst99 provides the tabulation of the values 
while GALEV does not. Therefore, individual spectra 
of GALEV SSP models of different ages were extracted 
and then Ha equivalent widths were measured using 
IRAF /splot with the same procedure as the one applied 
to observational data. The result comparing the two 
different SSP models is presented in Figure [2j For both 
Starburst99 and GALEV we assume an instantaneous 
burst of a stellar population with Salpeter IMF at solar 
metallicity, using Geneva stellar library. From the plot it 
is evident that HaEW may differ by a factor two or more 
between the two models. We attribute this difference 
partly because of the resolution difference of the model 
SED. GALEV SEDs have resolution of 20 A at Ha, 
thus the neighboring nitrogen lines ([Nll]AA6548,6583) 
were not resolved thus contaminated the Ha in the EW 
measurement. However, this EW difference translates 
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Table 1 

Target SN sites IFU observations 



SN 


Type 


RA2000 


Dec2000 




Galaxy (iVg N ) 


d/Mpc a 


Obs. date b 


Inst. c 


Exposure 


Seeing 


1964L 


Ic 


11:52:49.09 


+44:07:45.4 


±2" 


NGC 3938 (3) 


17.4 


2011 Mar 13 


S 


1800 


s 


x2 


0.8" 


19941 


Ic 


13:29:54.07 


+47:11:30.5 


±1" 


NGC 5194 (4) 


7.9 


2011 Mar 10 


s 


1800 


s 


x2 


1.1" 


2000ew 


Ic 


11:40:58.52 


+11:27:55.9 


±0.7" 


NGC 3810 (2) 


17.8 


2011 Mar 15 


s 


1800 


s 


x2 


0.9" 


2004gt 


Ic 


12:01:50.37 


-18:52:12.7 


±0.005" 


NGC 4038 (4) 


20.9 


2011 Mar 10 


s 


1800 


s 


x2 


1.2" 


2007gr 


Ic 


02:43:27.98 


+37:20:44.7 


±0.02" 


NGC 1058 (3) 


7.7 


2011 Sep 29 


G 


1800 


s 


x4 


0.5" 


2009em 


Ic 


00:34:44.53 


-08:23:57.6 


±1" 


NGC 157 (1) 


19.5 


2011 Sep 29 


G 


1800 


s 


x2 


0.6" 


1983N 


lb 


13:36:51.23 


-29:54:01.7 


±0.6" 


NGC 5236 (5) 


6.9 


2011 Mar 11 


S 


1800 


s 


xl 


1.5" 


1984L 


lb 


02:35:30.54 


-07:09:30.3 


±0.3" 


NGC 991 (1) 


18.8 


2011 Sep 29 


G 


2700 


s 


x2 


0.6" 


1999ec 


lb 


06:16:16.16 


-21:22:09.8 


±0.2" 


NGC 2207 (4) 


31.3 


2011 Mar 10 d 


S 


1800 


s 


x2 


1.1" 


2008bo 


lb 


18:19:54.41 


+74:34:21.0 


±0.2" 


NGC 6643 (2) 


22.1 


2010 Aug 1 


S 


1800 


s 


xl 


0.8" 


2009jf 


lb 


23:04:52.98 


+12:19:59.5 


±0.06" 


NGC 7479 (2) 


33.9 


2011 Sep 29 


G 


1800 


s 


x4 


0.5" 



a Mean redshift-independent distance from NED (http://ned.ipac.caltech.edu). 

b Hawaiian Standard Time (UTC 10). 

c Instrument used: S = SNIFS, G = GMOS. 

d Great East Japan Earthquake occurred during the exposure of this object; observation was conducted remotely from Tokyo but fortunately 
the earthquake did not introduce any effect to the observation and data. 
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Figure 2. Comparison of Ha equivalent width evolution from 
different SSP models, Starburst99 and GALEV. 

into only a small difference in age, in the order of 
1-2 Myr. This, eventually, will lead into difference of 
typically only about 20-30 % or lower in star initial 
mass, depending on the region of the age probed by 
Ha equivalent width. Therefore, we conclude that the 
selection of SSP model does not introduce significant 
sytematics to the final result of progenitor star initial 
mass. In this case, Starburst99 was preferred over 
GALEV since it has better temporal resolution. 



3. THE EXPLOSION SITES 
3.1. SNe Ic sites 

1. SN 1964L site 

SN 1964L exploded in NGC 3938, a spiral galaxy host- 
ing two other SNe (type II SNe 1961U and 2005ay - see 
Paper II). This SN is poorly studied and no reference 
discussing the nature of the progenitor of this particu- 
lar ob ject was found in the literature. iBlavlock et al.l 
(|20QQh recovered the photographic spectra of the SN and 
classified it as type-Ic. We found that the cited posi- 
tion of SN 1964L is within 3 arcsec of a bright cluster 
in the host galaxy. Based on the measurements of the 
Palomar survey plates , the position was shown to be ac- 
curate to within 2" bv lPorterl (j!993f ). SNIFS data shows 
that this source exhibits a blue continuum with strong 



emission lines, typical of a young stellar population. Fig- 
ure [3] shows the appearances of the explosion site on the 
IFU focal plane and the extracted spectrum of the host 
cluster. The measured metallicity is 1.09 solar. Mea- 
surement of Ha equivalent width shows that the cluster 
is very young, with derived age of 3.3 Myr. The corre- 
sponding stellar initial mass of that age at solar metal- 
licity is about 120 M . 

2. SN 19941 site 

SN 19941 is one of the best-studied type-Ic SNe due to its 
location in M51, a very well-observed nearby galaxy. It 
has often been referred as the 'standard' SN Ic although 
it al so exhibited unu s ual b ehaviour (e.g. iSauer et al.l 
I2006D . iNomoto etldl (|1994fl suggests that the progen- 
itor of this SN is a 4.0 M helium star in a binary sys- 
tem rather than a single WR star, which later evolved 
into a C+O star and eventually produced a SN Ic ex- 
plosion. The main sequence initial mass of this star was 
estimated as 15 M . The p osition of t his SN is accurate 
within 1", as estimated bv Ivan Dvkl (|1992l ) for SNe of 
that er a, and later confirmed b y the late-time observa- 
tions of iClocchiatti et all ([20081 ). 

Our SNIFS result shows that the explosion site is 
rather diffuse, the host cluster does not stand out promi- 
nently over the stellar background (Figure |4j). The ex- 
tracted spectrum shows faint continuum and a weak Ha 
line. Metallicity was derived to be 0.83 solar from N2 - 
this SN is the only SN Ic without an 03N2 determina- 
tion in our sample since the lines H/3 and [OIII]A5007 
are not detected. The Ha equivalent width corresponds 
to an age of 11.0 Myr. This age at solar metallicity cor- 
responds to the lifetime of a 17.9 M star, which is too 
low for a single WR progenitor but quite consistent with 
the binary model of INomoto et al.l (|1994l ). This is possi- 
bly one of the first strong pieces of evidence of a binary 
progenitor channel in SN Ic production. 

3. SN 2000ew site 

SN 2000ew exploded in NGC 3810 and is the first 
SN Ic with detected emission of near-infrared carbon 
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Figure 3. Left panel: SN 1964L site reconstructed IFU FoV. SN position within 1 arcsec error radius is indicated by a circle. Approximate 
linear scale corresponding to 2 arcsec is also indicated; this scale is calculated from the host galaxy distance only thus does not take into 
account the projection effect and host galaxy inclination. "SC" indicates the host star cluster. Right panel: extracted spectrum of the star 
cluster. 
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Figure 4. IFU FoV and extracted host cluster spectrum for the SN 19941 site. 
Figure annotations are the same as in Figure [3] 



The red end of the spectrum is truncated due to DAR. 



monoxide (CO ) in th e spectrum (Ge rardv et al.ll2002l ). 
I Van Dvk et al.l ([2003) searched the progenitor star in the 
pre-explosion images but found none within their detec- 
tion limit. Comparing with HST images taken while the 
SN is still visible, the derived SN position is accurate 
to within 0.7". However, they discovered that luminous 
blue-yellow stars are present in the explosion site, within 
1 arcsec of the SN position. The color-magnitude dia- 
gram suggests that the environment i s very young, < 6 
Myr old. Similarly, iMaund fc Smarttl (|2005D also failed 
to detect the progenitor star in pre-explosion images and 
derived an age estimate of the environment within 6 
arcsec from the SN of ^ 7 Myr, a ssuming twice solar 
metallicity. IMaund fc Smarttl (j2QQ5h estimated progeni- 
tor mass between 12 and 40 M from the environment 
age. 

Our SNIFS result shows that the cluster coincident 
with SN position, indicated as SC-B in Figure [5j is in- 
deed very young wi th age of 5.8 My r. This is consistent 
with the estimate of lVan Dyk et al.l (|2003f ). The brighter 
neighbouring cluster, SC-A, is slightly older, 6.3 Myr. 
Both clusters have similar metallicities of 1.15 and 1.05 
solar, respectively. This on-site metallicity measurement 
is definitely more reliable compared to the twice-solar 
metallicity assumed in the previous studies. The corre- 
sponding main sequence initial mass for a 5.8 Myr stel- 
lar population at solar metallicity is 33.9 M , while it is 
29.5 M for the 6.3 Myr population. 

I SN 2004gt site 



iGal-Yam et~aH ((2005), and also IMaund et al.l ([20051 ) re- 
ported their efforts in recovering the pre-explosion pro- 
genitor of SN 2004gt in the Antennae Galaxy, NGC 4038. 
However, the detection results were negative and only up- 
per luminosity limits could be derived for the progenitor 
star (positional uncertainty 5 mas). Both investigations 



could not rule out a massive Wolf-Rayet star progenitor, 
and also allow binary interaction scenario for producing 
the stripped-envelope progenitor star. The positional un- 
certainty of the SN in the HST images was estimated to 
be around 5 mas. 

Our result with SNIFS shows that SN 2004gt exploded 
in a large st ar cluster complex (Fig ure [6]), identified 
as Knot S by IWhitmore et al.l ([201(1 ). We derived the 
metallicity of Knot S as 12+log(0/H) = 8.71, or 1.12 so- 
lar, with age of 5.78 Myr. Th is result is consis t ent w ith 
Knot S mean age estimate by IWhitmore et al.l (|2010l ) of 
5.75 Myr. The age of 5.78 Myr corresponds to the life- 
time of a star with initial mass of 33.7 M at solar metal- 
licity. We also detect WR star feature in the spectra at 
4650 A. A metallicity determination of the explosion site 
bv lModiaz et~aTI (j20lTh gave 12+log(0/H) = 8.70±Hi> 
using the PP04 03N2 scale. Our 03N2 determination 
yields a higher value, 12+log(0/H) = 8.78. The dif- 
ference may have come from the different regions being 
observed, due to the different spectroscopy techniques 
(slit vs. IFU). However, if the uncertainty in the PP04 
03N2 calibration of 0.14 dex is considered as the errors 
for both determinations, this difference is statistically in- 
significant, corresponding to only 0.4a. 

5. SN 2001 r gr site 

We used GMOS to observe the explosion site of 
SN 2007gr which e xploded in the galaxy NGC 1058. 
iHunter et al.l (j2009[ ) reported an extensive observation 
of the SN at optical and infrared wavelengths, covering 
more than 400 days. They showed that the photometric 
evolution of SN 2007gr is similar to SN 2002ap, an ener- 
getic Ic hypernova, but the spectra of t he two SNe show 
marked differences. iTanaka et al.l (j2008[ ) performed spec- 
tropolarimetric observation of the SN and suggested that 
a bipolar explosion viewed slightly off-axis may best rep- 
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Figure 5. IFU FoV and extracted cluster spectra for the SN 2000ew site. Figure annotations are the same as in Figure [3] 
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Figure 6. IFU FoV and extracted host cluster spectrum for the SN 2004gt site. Figure annotations are the same as in Figure 03 



resent the explosion model of SN 2007gr. A study u sing 
pre-explosion imaging data bv lCrockett et al.l (|2QQ8al) re- 
ported that the SN exploded very near (~ 7 pc) a point- 
like source in HST images, which they proposed as the 
compact host star cluster of the SN progenitor star. The 
positional uncertainty of the SN is about 0.02" . Their de- 
termination of the compact cluster age using broadband 
photometry gives two solutions: 7 and 20-30 Myr assum- 
ing solar metallicity, corresponding to turnoff masses of 
28 and 12-9 M , respectively. 

GMOS observation of the explosion site was seeing- 
limited at 0.5". Despite the excellent seeing condition, 
the compact clusters and bright individual stars compris- 
ing the OB association resolved in HST images could not 
be resolved and they appear collectively as bright knots 
in the IFU reconstructed images. We extracted the spec- 
tra of each of these five sources, designated SC-A through 
SC-E, including the central association where SN 2007gr 
exploded (SC-A). We derived the metallicity of SC-A as 
1.12 solar, with age of 7.8 Myr, qui t e cons istent with 
the 7 Myr solution of lCrockett et al.l (|2008at ). This age 
corresponds to the lifetime of a 24.4 M star at solar 
metallicity. iModjaz et al.l (|2Qllh measured the metallic- 
ity at the explosion site to be 12 + log(0/H) = 8.64±q 09 
in PP04 03N2 scale. Our 03N2 determin ation yields 12 
+ log( 0/H) = 8.69, which is different with lModiaz et al.l 
(j2011f Ts result only by 0.25<j - this is not significantly 
different. 

Compared to other surrounding knots, SC-A appears 
to be older by about 1 Myr. All other knots SC-B 
through SC-E have ages between 6.4 and 6.8 Myr, which 
may suggest triggered star formation propagating inside- 
out from SC-A. If star formation in the outer clusters 
was triggered from SC-A, then the 1 Myr time difference 
and ~ 50 pc projected distance would correspond to a 
projected propagation velocity of ~ 50 km/s. The Ka 
map of the region shows shell-like structure surrounding 
SC-A and coincident with the other four knots, leaving 
a depression in Ha intensity at the position of SC-A - 
further reinforcing this scenario. In any case, the ages 



of the other four knots corresponds very well with the 
lifetimes of massive stars, around 28-29 M . Interest- 
ingly, SC-E shows relatively lower, sub-solar metallicity 
(0.85 Z ) compared to other knots present at the site, 
which are super-solar (> 1.1 Z ). 

6. SN 2009em site 

The Ic SN 2 009em is not well studied, discovered by 
iMonardl (|2009l) about 2 months after maxim um light as 
confirmed bv iNavasardvan fc Benettil (j2009f ) and also by 
iFolatelli fc Morrelll (j2QQ9h . The spectrum matches nor- 
mal type Ic events. There is no other information on 
this object available in the literature. We assign posi- 
tional uncertainty of the SN as 1", considering that it 
is a modern SN. The host galaxy NGC 157 is a grand- 
design Sc spiral, harbouring no other SN to date except 
2009em. 

With GMOS we observed the explosion site and de- 
tect two clusters. The more prominent cluster, SC-A, 
exhibits sub-solar metallicity of 0.76 Z , with age of 6.8 
Myr. The neighbouring fainter cluster SC-B is barely 
detected in continuum images but stands out in H<x It 
has nearly twice the metallicity of SC-A, 1.45 solar. The 
age of this cluster was derived as 6.3 Myr. The signifi- 
cant difference of the two clusters' metallicities illustrates 
the importance of measuring the metallicity of explosion 
sites with fine spatial resolution. The SN itself exploded 
at about the same distance from both clusters, but we as- 
signed SC-A as the host cluster since it is the brighter one 
and thus more likely to host more progenitor candidates. 
The ages of the clusters correspond to the lifetimes of 
27.7 and 29.5 M stars, for SC-A and SC-B respectively. 



3.2. SNe lb sites 



1. SN 1983N site 



SN 1983N exploded in a nearby large spiral galaxy, M83 
(NGC 5236). This SN was observed extensively (e.g. 
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Figure 7. IFU FoV and extracted cluster spectra for the SN 2007gr site. The continuum image of the GMOS FoV was made by collapsing 
the datacube in wavelength direction, approximately covering the V-band. Figure annotations are the same as in Figure [3] 
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Figure 8. IFU FoV and extracted cluster spectra for the SN 2009em site. Figure annotations are the same as in Figure [3] 
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Results for SN Ic sites 



SN site 


Object 


Offset/" (pc) a 


(03N2) 


(N2) 


12+log(0/H) 


Z/Z© 


HaEW/A 


Age/ My r 


M/M 


1964L 


SC* 


2.6 (220) 


8.83 


8.57 


8.70 


1.09tg-J? 


939.10 ±84.52 


o 9 o + 0.04 
O - ZO -0.04 


120^° 


19941 


sc* 


2.6 (110) 




8.58 


8.58 


0.83tg;£ 


13.34 ± 1.87 


11 o +0 - 25 


179-0.4 

1 ' -^+0.8 


2000ew 


SC-A 


2.3 (220) 


8.77 


8.58 


8.68 


l.UO_ 22 


68.18 ±8.18 


D - ol -0.04 


29.5;°;? 


2000ew 


SC-B* 


1.2 (110) 


8.85 


8.59 


8.72 


i - iO -0.30 


221.10 ±35.38 


^ 7C-+0.09 
°* /O -0.09 


33.9;°;* 


2004gt 


SC* 


1.3 (140) 


8.78 


8.64 


8.71 


19 +0.20 
1 - iZ -0.17 


209.60 ± 10.48 


c 7 o+0.03 
°- /o -0.03 


00 17 — 0.3 
°°- '+0.2 


2007gr 


SC-A* 


0.0 (0) 


8.69 


8.72 


8.71 


, -, 9 +0.03 
± - ±z -0.05 


15.57 ±0.96 


7 84+ - 75 


z ^'^+0.8 


2007gr 


SC-B 


1.4 (52) 


8.96 


8.63 


8.80 


1 oc:+0.64 
i - OO -0.42 


33.03 ±3.30 


6 66 +0 - 06 

°- DD -0.05 


OQ Q — 0.2 


2007gr 


SC-C 


1.1 (41) 


8.82 


8.68 


8.75 


1 20 +0 - 25 

±.^u_ 20 


59.62 ±2.47 




29.3;°;J 


2007gr 


SC-D 


1.6 (59) 


8.99 


8.58 


8.79 


1 op:+0.79 
i - OO -0.52 


50.26 ±5.03 


6 42+ - 06 

D -^ Z -0.04 




2007gr 


SC-E 


1.0 (37) 


8.65 


8.54 


8.60 


85 -0 - 14 


27.95 ± 1.39 


fi 7 o+0.08 


97 q 0-3 

z ' '^+0.2 


2009em 


SC-A* 


1.1 (100) 


8.57 


8.51 


8.54 


76+ - 05 

U - /D -0.05 


26.88 ±2.69 


fi o 4 +0.30 




2009em 


SC-B 


1.0 (94) 


8.85 


8.78 


8.82 


1 45+ - 10 


68.76 ±5.09 


D - O1 -0.03 


29.5-° l 



Offset between SN and approximate cluster center. 
SN parent cluster. 

Gaskell et all 119861 : IWeiler et all 119861 : iClocchiatti et all 



1996) and has been suggested to be the prototypical typ e 
lb SNe along with SN 1984L (jPorter fc Filippenkol[T987l ) . 
iSramek et al.l ([19841 ) reported t he first radio detection o f 
a type I SN of this SN 1983N. IClocchiatti etap (fl996h 
reported that their astrometry of the SN position is ac- 
curate to within ~ 0.6". 

Using SNIFS we discovered that the explosion site of 
the SN is quite complex, with three detected objects 
within the IFU FoV. In the continuum a prominent clus- 



ter was detected, but this cluster is overshadowed in Ha 
by two bright nearby Hi I regions (see Figure [9j). We 
extracted the spectra of SC-A cluster and the HI I re- 
gion north-east of it (SC-B), but failed to do so for the 
HI I region south-east of the cluster since it is situated at 
the edge of the FoV. It is likely that this third cluster is 
similar in age with SC-B considering its detectability in 
Ha and nondetection in continuum images. The cluster 
SC-A was found to be 7.2 Myr old with 1.02 times solar 
metallicity while the HI I region SC-B has age of 2.6 Myr 
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with 0.89 solar metallicity. The age of SC-A corresponds 
to the lifetime of a 26.4 M star while the SC-B age is 
considered too young to even produce a supernova. 

2. SN 1984L site 

SN 1984L explode d in NGC 991, a low surf ace bright- 
ness spiral galaxy. Schlegel & Kir shnerl (jl989l ) presented 
late-time photometry and spectroscopy of the SN, and 
sugges ted that it originat e d from a roughly 20 M star. 
Later, iSwartz fc Wheelerl (jl991l ) derived an eject a mass 
of more than 10 M©, which they infe rred to be the result 
of an explosion of a > 30 M star. iBaron et al.l ([1993) 
put forward the problem in reconciling the early-time 
light curve with the spectra of the SN, for which they pro- 
posed a solution by suggesting the scenario in which ex- 
plosion of a very massive star (~ 50 M of helium with 2 
x 10 52 erg kinetic energy) produced a massive black hole. 
The coordinates given by lSchlegel fc Kirshn er (1989 ) dif- 
fer by about 0.3" compared to those of IWeiler et al.l 
(1986). This value is used as the positional uncertainty 
of SN 1984L. 

Our GMOS observation of the explosion site shows 
that the host star cluster of SN 1984L does not exhibit 
strong emission lines in its spectrum, indicating interme- 
diate age. We measured the metallicity to be 0.45 solar, 
the lowest in our sample. Ha emission equivalent width 
is small, giving an age estimate of 18.0 Myr, which cor- 
responds to the lifetime of a 13.5 M star. This result is 
inconsistent with previous studies suggesting a high-mass 
(> 20 M ) progenitor, but intriguing since the derived 
mass of 13.5 M lies well below the currently accepted 
WR star mass limit and within the range of SN II-P 
progenitors. With such low metallicity, it is even more 
difficult to remove the hydrogen envelope of a single star 
via stellar wind mechanism, and as a result the minimum 
mass of a star to become a WR star increases. Thus, this 
is a strong indication that the SN 1984L progenitor may 
have lost its envelope in a close binary system. 

3. SN 1999ec site 

The lb SN 1999ec exploded in the interacting spiral 
galaxy, NGC 2207. Currently there have been four 
recorded SNe that have occurred in this galax y (SNe 
1975A, 1999ec, 2003H, 2013ai). IVan Dvk et "all ([200l ) 
searched for the progenitor star in pre-explosion images 
but found no convincing candidate. Their HST images 
show that the SN occurred in an environment of blue, 
apparently young stars and clusters. The SN position 
uncertainty was ±0.2". 

With SNIFS we detect two large clusters in the SN 
vicinity (Figure [TT]). Both clusters have similar spectra 
and are of similar metallicity, with SC-A being 0.71 solar 
value and SC-B of 0.74 solar. Also in derived age, both 
show similarity, with SC-A being 5.3 Myr and SC-B be- 
ing 5.1 Myr. These age estimates at solar metallicity 
correspond to the lifetimes of 38.0 and 41.6 M stars, 
respectively. The derived initial mass is high, even com- 
pared to the majority of SN Ic progenitors. 

I SN 200 8bo site 



SN 2008bo exploded in a spiral ho st, NGC 6643 . This 
SN has been detected in X-ray ( Immler" et al.l 2008; 
llmmler fc Brownll2008f ) and radio (Stockdale e tld1l2008[ ) 
wavelengths . The Asiago database as signed the type lb 
for this SN; iNavasardyan et aD (j2QQ8h reported that the 
spectrum shows the characteristics of a stripped-envelope 
event but with deep Ha absorption. This resembles the 
spectrum of SN 2008ax, a type lib SN, which shows spec- 
tral evolution from type II to lb. There is no other useful 
discussion about the physical natu re of the object in the 
literature. iStockdale et al.l (j2QQ8h reported that the ra- 
dio position of this SN agrees well (~ 0.2") with the 
optical position. Thus, we assign 0.2" as the positional 
uncertainty of this SN. 

We observed the explosion site using SNIFS and ex- 
tracted the spectra of two clusters apparent at the 
site (Figure [32]). A prominent cluster, SC-A, appears 
in all wavelengths while in Ha we detect a faint sec- 
ondary cluster, SC-B, and some other fainter structures 
in the FoV. We performed additional imaging observa- 
tion of the SN 2008bo s ite using OPTIC camera at UH88 
(jHowell fc Tonrvll2003f ) in 2011. However, the surround- 
ing fainter structures are not visible in our broadband 
images. From SNIFS data SC-A was found to be 13.5 
Myr old with 1.20 solar metallicity. We only managed 
to extract a part of the spectrum of SC-B but neverthe- 
less recovered the Ha region and derived an age via its 
equivalent width, which is 6.4 Myr. The metallicity of 
SC-B is derived as 0.85 solar. The 13.5 Myr age corre- 
sponds to a star with initial mass of 14.9 M©, while 6.4 
Myr corresponds to 29.3 M . The 14.9 M derived mass 
is too small for a star to produce a WR star that may 
later explode as a SN lb, thus it is likely that SN 2008bo 
also may have been produced by a progenitor in binary 
system. 

5. SN 2009jf site 

SN 2009jf exploded in one of the two prominent spi- 
ral arms of the barred spiral galaxy NGC 7479. It was 
an energetic slow -evolving SN, indicating massive ejecta. 
iSahu et all (|201lD presented 9-months of optical photom- 
etry and spectroscopy monitoring of this SN and sug- 
gested the main sequence mass of the progenitor to be 
> 20-25 Ma) estimated from the ejecta mass. Further, 
IValenti et al I (|2011f ) also reported their 1-year monitor- 
ing of SN 2009jf and similarly concluded that the pro- 
genitor mass range is within 25-30 M . Their positional 
uncertainty of the SN is in the order of 0.06". 

We observed the explosion site of SN 2009jf using 
GMOS, two years after the discovery. The contin- 
uum IFU reconstructed image shows only one prominent 
source in the fie ld (Figure [T3|). whi ch are the blended 
sources A+B in IValenti et al.l (j2011f ). We assigned the 
designation SC-A star cluster for this blended source. 
Making use of the capability of IFU we reconstructed 
Ha image of the field, and found that clumpy shell-like 
structure of ionized gas is present around SC-A. We as- 
signed identification SC-B, SC-C, and SC-D for these 
clumps. We found that the metallicity of this complex is 
below solar, with SC-A being 0.78 solar value and SC-B, 
C, D within 0.63-0.81 solar. The N2 metallicity of SC- 
A in 12+log(0/H) equals 8.55 dex, agree s exactly with 
the explosion site metallicity estimate of San ders et al.l 
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Figure 9. IFU FoV and extracted cluster spectra for the SN 1983N site. Figure annotations are the same as in Figure [3] The blue part 
of SC-B spectrum could not be extracted since it fell outside the FoV due to differential atmospheric refraction. 
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Figure 10. IFU FoV and extracted host cluster spectrum for the SN 1984L site. Figure annotations are the same as in Figure [3] 
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Figure 11. IFU FoV and extracted cluster spectra for the SN 1999ec site. Figure annotations are the same as in Figure [3] 
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Figure 12. IFU FoV and extracted cluster spectra for the SN 2008bo site. Figure annotations are the same as in Figure [3] Only a part 
of SC-B spectrum could be recovered due to DAR effects. 



(f2012h . using slit spectroscopy of 0.7" slit width. How- 
ever, they do not mention about the aperture size for 
spectru m extraction n o r the seeing size during the obser- 
vation. S anders et al.l (|2Q12h determined the age of the 
explosion site from H/3 equivalent width to be 9.0 Myr. 
The age of SC-A from our determination is quite old, 
18.2 Myr which corresponds to the lifetime of a 12.4 M 
star. The surrounding clusters are significantly younger, 
between 6.0 and 6.4 Myr. These morphological char- 
acteristics are similar with what we found at the SN 
2007gr explosion site, where the central host cluster in 
the complex is older than the surroundings thus suggest- 
ing inside-out star formation activity. T his formation 
scenario has been observed elsewhere (e.g. Ad amo et al.l 



l2012f h thus suggesting that triggered star formation may 
not be uncommon in star formation complexes. Assum- 
ing the outer clusters formation was triggered from SC-A, 
the projected outward velocity of star formation propa- 
gation in this complex is ~ 300 pc per 12 Myr, or around 
25 km/s. This is about half the projected velocity we in- 
ferred for the SN 2007gr site. 

If the progenitor of SN 2009jf indeed originated in SC- 
A, this would be a third example after SN 1984L and SN 
2008bo of low-mass SN progenitor produced by binary 
channel. In contrast, the surrounding environment age 
is consistent with the lifetime of stars with mass around 
30 Mq. We note that the progenitor mass inferred from 
SN properties is more consistent with the outer envi- 
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ronment age rather th an the stellar assoc iation closest 
with the SN position. IValenti et al.l (j2011f ) found that 
the color of their source A is consistent with a stellar 
population containing stars with maximum stellar mass 
of 8-25 M , while source B is redder thus presumably 
older. This is not inconsistent with our findings which 
yield progenitor mass of about 12 M . However, it is 
possible that the star formation history in the immedi- 
ate vicinity of SN 2009jf is not instantaneous, and the 
progenitor star was born from a younger burst compared 
to the dominant older surrounding population. 

4. DISCUSSIONS 

It is very important to determine the metallicity and 
age of the explosion sites of SNe Ib/c to characterize their 
progenitor stars. Metallicity-driven wind is one widely 
accepted scenario to remove the star's hydrogen envelope 
and produce a stripped SN. Our results show that even at 
the SN explosion site metallicity variations from place to 
place may be present, thus emphasizing the importance 
of measuring local metallicity directly with the highest 
spatial resolution possible and discouraging metallicity 
determinations using proxies. The results are presented 
in Tables [3] and El 

Averaging metallicities of all the 12 objects present at 
each explosion sites, we found mean metallicity for the 
sample of SN Ic sites as 1.11 ± 0.22 (RMS; standard er- 
ror of the mean/SEM = a/y/N = ±0.06) Z©, while this 
value is reduced to 1.01 ±0.17 (SEM = ±0.08) Z if only 
the SN parent clusters are considered. For the case of SN 
lb sites the mean metallicity value is 0.79 ± 0.20 (SEM 
= ±0.06) Z , and 0.83 ± 0.29 (SEM = ±0.13) Z for 
SN lb host clusters only. Notwithstanding the somewhat 
small number of samples (6 SN Ic sites, 12 clusters; 5 SN 
lb sites, 11 clusters), this immediately suggests that SNe 
Ic are produced in higher metallicity environments than 
SNe lb. The statistical significance of the difference be- 
tween SN Ic and lb host clusters metallicities determined 
by the t-test is 1.2<r, while it is 3.7a if we consider the 
si te metallicities . 

iModjaz et all ([201 lh argued that SNe Ic and lb ex- 
plode in different metallicity environments with a differ- 
ence of 0.2 dex in 12±log(0/H), with total number of 
sample of N = 35. Their Kolmogorov-Smirnov (KS) test 
p-value of the null hypothesis of both SN types came 
from the same population yields 1%. Converting back 
our solar-unit metallicity values of the host clusters into 
12+log(0/H) scale, we found 0.08 dex difference between 
SN Ic and lb explosion sites. The KS test p-value for our 
two samples is 59%, which could not rule out the possi- 
bility that both samples are from the same population - 
however this might not be very meaningful considering 
the s mall number stat i stics. For comparis on, other stud- 
ies by iLeloudas et al.l ([201 U N = 20) and lSanders et al.l 
(2012, N = 33) also suggest higher meta llicity for Ic (al- 
beit not statistically significant), while I Anderson et al.l 
(|2010l ) found no difference be tween lb and Ic meta llici- 
ties (N = 24). The sample of iSanders et al.l ([2012( 1 was 
drawn from untargeted SN surveys, thus is likely to be 
less biased towards massive, high- metallicity host galax- 
ies compared to a sample drawn from targeted SN sur- 
veys. Our result (N = 11) of 1.2a significance suggests 
that on average SNe Ic are the more metal rich compared 



to SNe lb, but not statistically si gnificant. There f ore, it 
is more closer to the conclusion of ILeloudas et al.l ()2Qlll ) 
who also found suggestion that SNe Ic are more metal 
rich compared to SNe lb, similarly by 0.08 dex. 

We note one interesting fact that for SN Ic sites the 
metallicity derived from 03N2 index is almost always (10 
out of 11) higher than from N2 index (3.7a significance), 
while it is the opposite for SN lb sites (2.7a significance). 
Consistently for both SN Ic and lb sites the measured 
metallicity is always higher for SN Ic sites in either 03N2 
or N2 scales. The high Q3 N2 metallicity may be in- 
dicative of N/O enhancement ( Per ez- Mont ero fc ContinU 
2009). This behavior has been observ ed in galaxies wit h 
high 03N2 metallicity as reported by lBerg et all ([201 lh . 
and interpreted as galaxies evolving past Wolf-Rayet 
galaxy phase as WR stars may elevate the observed N/O 
ratio. 

The age of the SN explosion site provides clue to the 
likely mass of the progenitor star. We averaged the age 
of all detected clusters in our SN sites and found that SN 
Ic sites have younger mean age than SN lb sites, 6.6 ± 1.7 
(SEM = ±0.5) Myr compared to 8.6 ±5.3 (SEM = ±1.6) 
Myr. This trend is also present when comparing the ages 
of the host clusters, with SN Ic hosts averaged 6.7 ± 2.6 
(SEM = ±1.1) Myr while for SN lb hosts the age is sig- 
nificantly higher, 12.4 ± 5.9 (SEM = ±2.6) Myr. This 
result shows that it is likely that SNe Ic resulted from 
higher-mass progenitors than SNe lb (2.0a significance). 
Considering the results for individual objects, we also 
suggest that some SN Ib/c progenitors may have origi- 
nated from close binary systems due to the derived old 
age of the host cluster. SN 19941 is one example for 
SN Ic with derived old parent cluster age, or low mass 
progenitor. The derived age of 10.9 Myr corresponds 
to the lifetime of a 17.9 M single star, which is not 
massive enough to evolve into a stripped-envelope WR 
star thought necessary to produce SN Ic. Therefore, one 
likely scenario is that this SN was the explosion of a star 
in a binary system, whose hydrogen envelope has been 
removed by binary interactions. For SNe lb in our sam- 
ple we found that this scenario may be applicable to SNe 
1984L, 2008bo, and 2009jf. It is also interesting to note 
that the binary scenario is significantly more prevalent 
in our SNe lb samples (3 out of 5 SNe) compared to SNe 
Ic (1 out of 6 SNe), but further investigation with more 
samples is necessary to confirm this indication firmly. 

It is interesting to note that type lib SNe, whose spec- 
tral evolution showed transition from type II to type lb, 
already show evidence th at points to a s t ar in a binary 
system as the progenitor. lAldering et al.l (j!994f ) showed 
that the SED of SN 1993J proge nitor could no t be fit 
with a single star SED; further, iMaund et all (j2004f ) 
reported the unambiguous detection of the companion 
star with post-explosion observation. Indication of bi- 
narity was also found in recent lib event s: SN 2008ax 
(iCrockett et al.l l2008bf ) and SN 2011dh ([Bersten et al.l 
l2012f ). Based on the similarities, it is possible that both 
type lib and lb SNe have similar binary progenitors. 

With the information on SN host cluster age and 
metallicity, we derive the maximum mass of the stars 
still present in the cluster. As the most massive stars 
dies first, any of the remaining stars would not exceed the 
mass of the SN progenitor exploded just recently. The 
age of the cluster is equal to the lifetime of the SN progen- 
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Figure 13. IFU FoV and extracted cluster spectra for the SN 2009jf site. Figure annotations are the same as in Figure [3] 



Table 3 

Results for SN lb sites 



SNsite Object Offset/" (pc) a (03N2) (N2) 12+log(0/H) Z/Z HaEW/A Age/Myr M/M 



1983N 


SC-A* 


2.6 (57) 




8.67 


8.67 


1.02 


1983N 


SC-B 


3.0 (66) 




8.61 


8.61 


0.89 


1984L 


SC* 


1.6 (150) 




8.31 


8.31 


0.45 


1999ec 


SC-A* 


0.8 (90) 


8.49 


8.53 


8.51 


0.71 


1999ec 


SC-B 


3.3 (370) 


8.51 


8.55 


8.53 


0.74 


2008bo 


SC-A* 


3.1 (330) 




8.74 


8.74 


1.20 


2008bo 


SC-B 


3.5 (380) 




8.59 


8.59 


0.85 


2009jf 


SC-A* 


0.0 (0) 




8.55 


8.55 


0.78 


2009jf 


SC-B 


1.4 (230) 


8.42 


8.53 


8.48 


0.66 


2009jf 


SC-C 


1.0 (160) 


8.39 


8.52 


8.46 


0.63 


2009jf 


SC-D 


1.2 (200) 


8.50 


8.63 


8.57 


0.81 



+ 0.53 
-0.34 
+ 0.46 
-0.30 
+ 0.13 
-0.16 
+ 0.30 
-0.30 
+ 0.04 
-0.02 
+ 0.62 
-0.41 
+ 0.44 
-0.29 
+ 0.39 
-0.27 
+ 0.08 
-0.08 
+ 0.09 
-0.09 
+ 0.12 
-0.12 



23.32 ±4.19 
1693 ± 728 
2.93 ±0.49 
384.50 ±46.14 
424.10 ±72.09 

3.13 ±0.53 
56.37 ±48.48 

1.38 ±0.12 
140.6 ±23.90 
88.71 ±9.76 
80.36 ±9.64 



7.22 
2.62^ 



+0.27 
-0.42 
) + 0.65 
-1.26 

18.012 2? 

5-34t° i 5 5 
5.06tr r ? 

iq k + 1.03 

lO.U_ 4Q 

6.37+ 5 - 92 
18.2+ 



-0.24 
>+0.77 
-0.52 
6-02 + 0.16 



-0.16 
. + 0.13 
-0.08 

6 2^+ - 04 
D - zo -0.05 



6.36^ 



26.4: 



13.5 
38.0 
41.6; 
14.9 



29.3 



+1 d 

- + 4 5 8 C 
+ 14. 
-0.1 

-13. 
+0.9 



-0.1 
*+0.f 

' -J-l.s 

-l.( 
+o.e 
-o.s 

+ 0.] 



12.4T 



32.2 
29.7 



a Offset between SN and approximate cluster center. 
SN parent cluster. 

itor, and since the lifetime of a single star is mainly gov- 
erned by its initial mass, the lifetime of the SN progeni- 
tor could be used to estimate its initial mass. With this 
method we determined the main-sequence mass of the SN 
progenitors in our sample by comparing to Padova stellar 
evolution models, adopting the metallicity of t he pa rent 
cluster. We used the models bvlBr essan et al. (1993|) for 
solar metallicity and iFagotto et al.l (|1994f ) for 0.4 solar 
metallicity. The resulting mass along with metallicity 
information for each SN progenitor is plotted on mass- 
metallicity space in Figure [HI Theoret ical predictions 
of SN progenitors bv lGeorgv et al.l (|2QQ9[ ) are overlaid to 
make the comparison readily available. It is immediately 
apparent that the theoretical models and observational 
data points do not agree very well - only one occurrence 
of observational data falls into the right theoretical pre- 
diction for each Ib/c type (SNe 1983N and 1964L). In 
spite of this, we argue that the relative comparison of 
the mass of progenitor star of different SN types (i.e. 
Ib versus Ic) should be robust regardless of the stellar 
evolution model used. We also compare our results of 
progenitor mass and local environment age with other 
observational results using different methods, summa- 
rized in Table |4j In general our results do not contra- 



2 
5 
i 




diet the findings of other studies, except probably for 
SN 1984L. This provides a strong confirmation of the 
validity of our method. However, we note that there is 
always probability that we are suffering from projection 
effects which may cause apparent association between 
the SN and an otherwise unrelated stellar population, 
resulting in a chance superposition. This issue is fur- 
ther addressed in Paper II, resulting in a 50% chance 
superposition estimate. Nevertheless our method offers 
the current best determination of SN progenitor age and 
metallicity from the SN local environment. Median dis- 
tance of objects in this study is ~15 Mpc (1" corre- 
sponds to 70 pc), or redshift 0.0035 (for Hq = 70 km s _1 
Mpc -1 ), the smallest compared to other similar studies. 
For example the medi an redshift (distance ) of objects in 
lAnderson etaLl (j2010D is 0.005 (21 MpcUModjaz et al.l 
([201 lh Qj)lf772 Mpc ^ fLeloudas et all ((2011|) 0.022 (94 
Mpc), [Sa nders et~aU (|20ll 0.036 (154 Mpc) (Table 8 of 
iSanders et al.ll2012l ). 

In this work we suggest that SNe Ic originate from 
more massive progenitors than SNe Ib. The progenitor 
mass of SN Ic averages 42.9 ± 38.2 M , while for SN 
Ib the average progenitor mass is 21.0 ± 11.0 M . The 
high 03N2 metallicity of SN Ic sites supports this con- 



SN site IFU spectroscopy - Ib/c progenitors 



13 



2.0 




99ec 



0.0 



~i — i — i — i — i — i — r 



SN Ic 



▲ IC 




64L 







20 40 60 SO 

M/Mq 



100 



120 



Figure 14. SN Ib/c progenitors overplotted on mass-metallicity space; Ic progenitors are indica ted with purpl e trian gles and lb progenitors 
are indicated with green diamonds. Theoretical predictions of progenitors of each SN type by Georgv et al. (2009) are drawn as colored 
regions. 

elusion that SN Ic progenitors have mass within the WR 
star range. We note that the large deviation in SN Ic 
average progenitor mass is caused by the very high de- 
termined mass of SN 1964L, ~ 120 M . If we omit this 
one SN the average progenitor mass for SN Ic is reduced 
to 27.5 ± 6.7 M - still higher than SN lb progenitor av- 
erage and consistent with WR star mass. The statistical 
significance of the lb and Ic progenitor mass difference 
is 1.3a, and the inclusion of SN 1964L will slightly re- 
duce it into 1.1a. Our results also show some examples 
where the progenitor may have originated from close bi- 
nary systems, considering the old age of the parent stel- 
lar population. This is consistent with the findings of 
iLeloudas et al.l (|2011f h which suggests that both single 
or binary progenitor scenarios are at work for SN Ib/c 



production. Recent finding by iSana et all (j2Q12h shows 
that binary interaction may affect 70% or more of the 
massive star population. 

Considering the two channels of high and low mass 
progenitors possible for SN Ib/c, it is interesting to com- 
pare it with the derived initial masses of the compact 
remnants in the Galaxy. Magnetars are believed to be 
formed by rapidly rotating massive stellar cores, and this 
mechanism would not work if the star enters the red 
supergiant phase where magnetic breaking between the 
stellar core and envelope may spin down the core rota- 
tion. Therefore, a very high mass star (> 40 M ) that 
can lose a significant portion of its hydrogen envelope 
prior to SN is required as the progenitor of magnetars. 
Several studies have confirmed this by studying the host 



14 



Kuncarayakti et al. 



Table 4 

Mass and age of SN Ib/c progenitors from this work, compared with other determinations 





SN 


Mass (this w.) 


Age (this w.) 


Mass (ref.) 


Age (ref.) 


Method (ref.) 


Reference 


Ic 


1964L 


120 M 


3.28 Myr 










Ic 


19941 


17.9 M 




15 M 




SN properties 


Nomoto et al. (1994) 


Ic 


2000ew 




5.75 Myr 




< 6 Myr 


environment CMD 


Van Dvk et al. (2003) 






33.9 M 




12-40 M 




nondetection 


Maund & Smartt (2005) 


Ic 


2004gt 




5.78 Myr 




5.75 Myr 


SC color 


Whitmore et al (2010) 






33.7 M 




> 40 M , or 




nondetection 


Maund et al. (2005), 








low-mass binary 






Gal- Yam et al. (2005) 


Ic 


2007gr 


24.4 M 


7.84 Myr 


28 or 12-9 M 


7 or 20-30 Myr 


SC color 


Crockett et al. (2008a) 


Ic 


2009em 


27.7 M 


6.84 Myr 










lb 


1983N 


26.4 M 


7.22 Myr 










lb 


1984L 


13.5 M 


18.0 Myr 


20 M 




SN properties 


SchleEel & Kirshner (1989) 



lb 1999ec 
lb 2008bo 
lb 2009jf 



38.0 M 
14.9 M 

12.4-34.7 M 
12.4-34.7 M 
12.4-34.7 M 



5.34 Myr 
13.5 Myr 
6-18 Myr 



> 30 M 

> 50 M 
My > -8.7 



>20-25 M 
25-30 M 
<8-25 M 



9.0 Myr 



SN properties 
SN properties 
nondetection 

site H/3EW 
SN properties 
SN properties 
SC color 



Swartz & Wheeler 
Baron et al. (1993) 
Van Dvk et al. (2003) 




Sand ers et al (2012) 
Sahu et aT720TlT 
Valenti et al. (2011) 
Valenti et al. (2011) 



Note. — CMD: color-magnitude diagram; SC: star cluster 

star clusters of the m agnetars (e.g. iMuno et al.l 12006), 
but lDavies et al.l (j2009[ ) showed that the magnet ar SGR 
1900+14 may have had a low mass progenitor of around 
17 M . Therefore, the SN that produced those magne- 
tars should also came from both massive and lower-mass 
stars, and the SN needs to be a stripped-envelope event 
to be able to produce rapid core rotation that eventually 
gave birth to the magnetar. In context with SN Ib/c 
progenitor initial mass derived from our study, the two 
results are in line. More samples of both Ib/c SN and 
magnetar progenitors would be beneficial to establish a 
firm connection between the two. 

We note that this study is based on the assumption 
that the SN progenitor had been a member of the appar- 
ent host star cluster present at the explosion site. While 
this assumption might not necessarily be true, we argue 
that it is more likely that the associated cluster is the 
real parent stellar population of the SN progenitor rather 
than the diffuse background stellar population. The host 
cluster has been revealed to have very young age, thus the 
IMF implies higher probability of having high- mass stars 
compared to the background stellar population which is 
significantly older than the few Myr age of the prominent 
cluster. We discuss the probability of contamination by 
older, invisible clusters in the field in Paper II, where 
the effect is more important for type II SN progenitors, 
which are presumably older than SN Ib/c progenitors. 

We also present our mass estimate if we only consider 
the high-association cases only, i.e. SNe within 150 pc 
of the host cluster center. This 150 pc is about the in- 
termediate size between classical and giant HI I regions 
(jCrow ther 2013), as most of our clusters are apparently 
of around this size. With this limitation, the mean pro- 
genitor mass of SNe Ic reduces to 27.5 ± 6.7 M , and SN 
lb progenitors to 22.6 ± 12.1 M . This would not change 
the aforementioned conclusion. 

5. SUMMARY 

We present the results of our investigation of nearby 
SNe Ib/c explosion sites using integral field spectroscopy 



(IFS) technique, performed using UH88/SNIFS and 
Gemini/GMOS. Taking advantage of IFS, we spatially 
identified star clusters at the explosion sites and ex- 
tracted their respective spectra for the purpose of analy- 
sis. From the spectra we derived the metallicities of the 
clusters in terms of 12+log(0/H) using the strong line 
method, and age using comparison of observed Ha emis- 
sion equivalent width with theoretical SSP models from 
Starburst99. 

The method employed in this study offers the best 
available determination of SN progenitor natal metallic- 
ity, probing a very localized region specifically the par- 
ent stellar population from which the SN progenitor was 
born. This minimizes the effect of contamination from 
other populations and definitely is more reliable com- 
pared to metallicity determinations using proxies. As 
metallicity is one important parameter in context of SNe 
Ib/c progenitor studies, our result provides an impor- 
tant insight in this field. We found indication that the 
metallicity of SNe Ic explosion sites is on average higher 
than SNe lb s ites, by ~ 0.1 d ex. T his r esult is in ac- 
cordan ce with iLeloudas et al.l (|2Qllh and San ders et al.l 
(j2Q12h who also found indication that SN Ic explosion 
sites are more metal-rich compared to S N lb sites, al- 
thoug h not statistically significant, and Mo djaz et al.l 
(j2011f ) who suggest the same conclusion with results they 
claim to be statistically significant. 

The age of the SN progenitor host cluster is consid- 
ered as the lifetime of the progenitor star. Comparing 
with Padova stellar evolution models, we are able to de- 
rive initial mass of SNe Ib/c progenitors from the ages of 
their parent star cluster. We found that SNe Ic explosion 
sites have younger age than SNe lb, indicating that they 
arise from higher mass progenitors than SNe lb. The 
derived progenitor mass from host cluster age is signifi- 
cantly higher for SN Ic, over 42 M versus 21 M for SN 
lb progenitors. We note that these derived average pro- 
genitor mass were derived from a sample that includes 
sub-WR mass progenitors (< 25 M ) which may be in- 
terpreted as binary progenitors. If we consider the sin- 



SN site IFU spectroscopy - Ib/c progenitors 



15 



gle progenitors only, SN Ic progenitors average 47.9 M 
while SN lb progenitor mass average is 32.2 M - still 
attesting that SN Ic progenitors are the more massive 
ones. Further, if only SNe showing high association with 
the parent clusters are considered, the average progen- 
itor mass estimates would reduce to 27.5 and 22.6, for 
SNe Ic and lb respectively. 

In general, our findings support the current picture of 
stripped envelope (Ib/c) SN progenitors: SNe Ic tend to 
have more massive progenitors and more metal-rich than 
SNe lb, and for both SN types interacting massive binary 
stars may comprise a fraction of the progenitors of the 
SNe. 
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